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Abstract 



We study the effect of the Higgs-exchange diagram for the lepton flavor violating 
muon-electron conversion process in nuclei in the supersymmetric seesaw model. The con- 
tribution is significant for a large value of tan f3 and a small value of a neutral heavy Higgs 
boson mass, in which case the ratio of the branching ratios of B(fiN — > eN)/B(fi — > ej) 
is enhanced. We also show that the target atom dependence of the conversion branching 
ratio provides information on the size of the Higgs exchange diagram. 
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Among various candidates for physics beyond the standard model, the supersymmetric 
(SUSY) extension is considered to be one of the most promising ones, providing us with 
a solution to the hierarchy problem. In addition, the gauge coupling constants measured 
precisely in the last decade show a remarkable agreement with the prediction of the SUSY 
grand unified theory (GUT) p. Although direct searches for SUSY particles are the most 
important, it is also interesting to see implications of such a theory in the low-energy 
phenomena. This will be important even after we discover SUSY particles at energy 
frontier experiments, because some of model parameters can be only accessible from low 
energy experiments. 

Lepton flavor violation (LFV) offers a possibility to explore SUSY models from low 
energy experiments In the standard model, the lepton number is conserved sepa- 
rately for each generation, so that LFV in charged lepton precesses is forbidden. On the 
other hand, recent developments in neutrino physics indicate that the lepton flavor is not 
conserved in the neutrino sector. However, the simplest model of a finite neutrino mass, 
namely the see-saw model, does not induce observable effects of LFV in charged lepton 
processes since only the mass terms of the neutrinos violate the lepton flavor conservation 
PJ. The situation is completely different in the context of SUSY. In SUSY models the 
LFV in muon and tau decays is considerably enhanced due to the existence of the scalar 
partner of leptons, and therefore the branching ratios of these processes can be close to the 
reaches of current or near future experiments in SUSY-GUT 0] as well as SUSY seesaw 
models UUU\- 

In the SUSY model, a new source of LFV appears in the off-diagonal components of 
the slepton mass matrices. In the case of the SUSY seesaw model, the Dirac Yukawa 
interactions of the neutrinos induce the off-diagonal components at the one-loop level, 
even if we assume the slepton mass matrix to be proportional to the unit matrix at the 
high-energy scale such as in the minimal supergravity scenario p]. This effect can be 
sizable, and the current upper bound of the branching ratio of the fi — > e-y decay already 
puts severe constraints on the SUSY parameters. 

Another important process is the fi — e conversion in nuclei. In the effective Lagrangian 
at the energy scale of the muon mass, this process can be induced by several four-fermion 
operators, in addition to the photon dipole-type operator, which is responsible for the 
fi — > e7 decay. If the latter is the only source of LFV, the branching ratio for the 
fi — e conversion is suppressed roughly by 0(a) compared to the branching ratio of the 
fi — > e7 decay. Even if this is the case, significance on new physics search from two 
processes can be similar since the experimental upper limit is lower for the /i — e conversion 
process. The current best experimental upper bounds for the branching fractions are 
B(fi -> ej) < 1.2 x 1(T U |H] and B(jiTi -> eTi) < 6.1 x 10~ 13 0, respectively. There are 
several planned experiments which are aiming at improving the bounds of the branching 
fractions for relevant processes by three or four orders of magnitudes fUJ El E3 • If other 
four-fermion interaction is sizable, or even dominant, the fi — e conversion branching ratio 
may not be suppressed by 0(a) relative to the /x — > branching ratio. For example, 
in R-parity violating SUSY models, the contribution from the scalar type four-fermion 
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interaction is shown to be important especially through the strange-quark and the 
bottom-quark couplings (Ej. 

Recently, the effect of the neutral Higgs exchange diagrams in the various flavor chang- 
ing neutral current (FCNC) processes is considered, and a possibility of large contributions 
is pointed out especially for large tan/? and small region in SUSY models In 
LFV processes, the new effect to the r — > 3/i and r — > fir] decay is studied [IBJ HZl UBJ- 
Since the Higgs- mediated FCNC does not contribute to the r — > /X7 decay, the ratio of 
B(t — > 3/jl)/B(t — > fij) is useful to reveal the existence of the effect. 

In this letter, we studied the effect of the Higgs-exchange diagram on the \i — e con- 
version process in the SUSY seesaw model. In contrast to the \i — > 3e decay, the Higgs- 
mediated contribution to the \i — e conversion process is not suppressed by the electron 
mass but only by the nucleon masses, because the Higgs-boson coupling to the nucleon is 
shown to be characterized by the nucleon mass using the conformal anomaly relation [19 . 
The most important contribution turns out to come from the exchange of the heavier scalar 
Higgs boson (H°) which couples to the strange quark scalar current in the nucleus. We 
found that the transition amplitude from this type of diagrams becomes fairly large com- 
pared to the photon-exchange diagram responsible for the \i —>■ ej decay in the large tan /3 
and the light H° region. Therefore, the ratio of B(fiN —>■ eN)/B(fi —>■ ej) is quite sensi- 
tive to the Higgs-exchange effect, just as B{r — > 3fi)/B(r — > fij) is important in r decays. 
Also, we show that it is possible to identify the Higgs-mediated LFV effect by looking at 
the target atom dependence of the branching ratio, e.g. B(fiPb — > ePb)/5(/xAl — > eAl). 

In the SUSY seesaw model, the off-diagonal components of the slepton mass matrix 
appear in the left-handed sleptons through the neutrino Yukawa interactions. We assume 
that the slepton mass matrix is proportional to the unit matrix at the GUT scale, and 
evaluate the effects of the neutrino Yukawa interaction to the slepton sector. The super- 
potential of the lepton sector is given by W = f e Ei ■ EfLi + f l u j H 2 ■ NfLj + (l/2)M l N N£ N[ , 
where Hi and H 2 are the doublet Higgs fields, L iy E?, and Nf are the superfields cor- 
responding to the left-handed leptons, right-handed leptons, and right-handed neutrinos 
of the i-th generation, respectively. The neutrino mass matrix is obtained by integrating 
out the heavy right-handed neutrinos as rrfj = (f^M^f^ijv 2 sin 2 (3/2 (2Uj, where v is 
the vacuum expectation value (VEV) of the Higgs field (v = 246 GeV) and the angle 
(3 is defined by tan (3 = (H®) / (H®) . To obtain the correct size of the neutrino masses, 
the right-handed neutrinos should be as heavy as 10 14 GeV for f v ~ 0(1). The Yukawa 
interactions represented by f v violate the lepton flavor conservation. This violation is 
imprinted to the slepton mass matrix in the low-energy Lagrangian. The renormaliza- 
tion group equation (RGE) running effect induces the off-diagonal components in the 
left-handed slepton mass matrix which are approximately given as follows: 

(Am?J,, ~ -^fprriKZ + |a | 2 ) log ^ , (1) 

where the SUSY breaking parameters mo and ao represent scalar masses and three point 
scalar interactions at the GUT scale, respectively. 
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LFV in the Higgs coupling originates from the non-holomorphic correction to the 
Yukawa interactions of the charged leptons ^B] • One-loop diagrams mediated by sleptons 
induce the following Yukawa interaction terms: 

£ = f&P^Hl + f% (4% + ef ] ) P L e 3 H Q 2 * + h.c. , (2) 
where Pl,r = (1 -F 7s)/2- The non-holomorphic interactions €\ and 6% are given by 1 



g 2 Y[ iM x 
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where g>y and g 2 are the gauge coupling constants, and Mi, M 2 , and /i are the gaugino and 
Higgsino mass parameters. The above formulas are based on the calculation of the effective 
Yukawa interaction in the SU(2)x, x U(l)y symmetric limit. The mass parameters tu\ 
and m~ are the slepton masses for the z-th generation. The functions 1% and I4 are 



defined by 

h(a,b, c) 



1 ab log(a/6) + be log(6/ c) + ca log(c/ a) 



(4tt) ; 



(a — b)(b — c)(c — a) 



(5) 
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+ 



(a — c) (6 — c) (d — c) (a — d) (b — d)(c — d) 



(6) 



Note that the parameters e\ and e~2 do not vanish even in the limit of large masses of SUSY 
particles. This is quite different from the photon-exchange diagrams of LFV, where the 
amplitude becomes small for large masses of internal SUSY particles. 

The Yukawa interactions in Eq.(J2J can be written in terms of the fields in the mass 
eigenstates. For the /1 — e transition, the Lagrangian is given by 



C 



(fiP L e) [cos(a - P)h° + sin(a - P)H° - iA°] + h.c. , 



(7) 



v cos 2 P 

where h° and H° are the scalar Higgs fields (m h o < m H o), and A is the pseudoscalar 
Higgs field. The LFV parameter k 2 \ is given by k 2 \ = 6^/(1 + tan/5) 2 . In the limit 

1 There are sign differences in theses equations compared to the results in Refs. |16l I17| . 
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where the masses of H° and A go to infinity, the LFV interaction of the lightest Higgs 
boson vanishes since the standard model does not have LFV. Therefore the contributions 
from H° and A are important in LFV processes for relatively small values of heavy Higgs 
boson masses. 

The n — e conversion process occurs through this interaction by the exchange of the 
Higgs boson with a nucleus. The effective four-fermion interactions are given by 



cff 



v cos 2 (3 



E 
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v sin (3 



cos(a — 0) cos a sin(a — (3) sin a 
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m 



H° 



(eP R fi)(qq) 
(8) 



The transition amplitude can be obtained by taking a matrix element. We evaluate the 
amplitude of the coherent conversion processes where the initial and final nuclei are in 
the ground state. Compared to incoherent transition processes, the coherent processes 
are expected to be enhanced by a factor of O(Z) where Z is the atomic number. In those 
processes, the matrix elements for the quark operators are obtained in the following way 
[2T] . The first step is to write down the effective Lagrangian in the nucleon level which is 
given by replacements of 



qq — > pp + G s nn 



and 



q^ 5 q -> G p y 'py 5 p + Up 'nj 5 n 



(9) 



where G's are coefficients which can be evaluated by taking matrix elements of quark 
operators by nucleon states. Then we can take the matrix elements by a specific nucleus. 
Since the initial and final states are the same, the elements (N\pp\N) and (N\nn\N) are 
nothing but the proton and the neutron densities in a nucleus in the non-relativistic limit 
of nucleons. In this limit, the other matrix elements (iV|p75p|iV) and (N\n , f^n\N) vanish. 
Therefore in the coherent \x — e conversion process, the dominant contribution comes from 
the exchange of H°, not A . 

As we can see in Eq.fjSJ), the LFV interactions contain factors of the muon and the 
quark masses because we need two chirality flips to form the scalar four-fermion operators. 
However, there can be also enhancement factors. One is the tan (3 enhancement in the 
H° vertices, with which the amplitude is proportional to tan 3 /3. The other factor of 
enhancement is due to the couplings of the Higgs boson to the nucleons represented by 



the G^' pS> and G K g' ll> factors. It is important to notice that the Higgs boson coupling to 
the nucleon is not suppressed by the up or down current quark masses because it can 
strongly couple to the gluons in the nucleon through the loop diagrams of the quarks [19 j . 
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Among the various quarks, the strange quark gives the dominant contribution in the large 
tan j3 region, since the down-type quarks have tan (3 enhancement in the Yukawa coupling 
constants. The values of the combinations of m s G^ and m s Gg' n ^ turn out to be much 
larger compared to the contribution from the down quark and that from bottom quark 
diagrams. The values are estimated to be rridGg /m p = 0.029, m d Gf n) /m n = 0.037, 
m s G { s' p) /m p = m s G ( s s ' n) /m n = 0.21, and m h Gf p) /m p = m b Gf n) /m n = 0.055 [22] • 

Once we write down the effective Lagrangian in the nucleon level, the estimation of the 
coherent conversion rate is straightforward [231 121] • The general interaction Lagrangian 
for the coherent conversion process is given by 

Ant = {m^A R Jia^P L eF^ + m^A L Jia^ P R eF^ + h.c.) 



Gp \ ^ 



g%ePRfJL + g l $eP Lf i)iH> 



+ 



~g^e^P L fx + g^e-fPsit) ^ + h.c. 



(10) 



where A's and g's are dimensionless coupling constants. The first two are the dipole 
operators which are the same operators for the fx — > decay. We also have scalar and 
vector type four-fermion operators. In the SUSY seesaw model, since the slepton mixing 
appears only in the left-handed sleptons, the parity is maximally violated, i.e. Al, g^ R s\ 
9rv^ ^ Ar ; 9ls ) i 9lv ■ All those operators are given through the one- loop diagrams 
mediated by the SUSY particles [U] • Besides the scalar operator from the Higgs exchange, 
only the dipole operator is important in the large tan j3 region because of the dependence 
of A R oc tan f3. 

With the above coefficients in the effective Lagrangian, the conversion rate is simply 
given by 



^conv — 2G F 



A* R D + 



111 



where D, V^, and are the overlap integrals among wave functions of 

the muon and the electron and the nucleon densities or the electric field in the nuclei. 
For example, in the aluminum nuclei, they are estimated to be D = 0.03Q2m 5 J 2 , = 
0.0155ttv /2 , = 0.0161t7v /2 , V® = 0.0167m^ /2 , and = 0.0173m^ /2 [21]. 

By comparing Eq.fJHJ) and Eq.Q with Eq. ljlUjl . we obtain the coupling constants g^s 
and g£gi and then we can derive the formula for the the conversion rate by Eq. ffTT]) . For 
example, in the aluminum and lead targets, the conversion branching ratios at large tan j3 
are approximately given by 



7 2 2 

B{fxM -> eAl) ~ 1.8 x 10~ 4 ■ ^W" 21 t an 6 /3 , (12) 
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and 



B(//Pb -> ePb) ~ 2.5 x 10 




tan b /3 , 



(13) 



v 4 m^ ix> ca pt 



respectively, where cu ca pt is the muon capture rate in the nuclei. The values are uj ca , pt = 
0.7054 x 10 6 s _1 and 13.45 x 10 6 s _1 in the aluminum and the lead nuclei, respectively [25 . 
If we take all the right-handed neutrino masses to be 10 14 GeV, the contribution of the 
Higgs exchange in Eq.(|12j) is roughly given by 



where Ms is defined by Ms = m = M\m with the universal scalar mass m and the 
gaugino mass M%/2 at the GUT scale. The above estimation shows that the Higgs exchange 
is important for tan/5 > 60 and m#o <C Ms region because of the different decoupling 
behavior. 

We numerically calculate the branching ratios in order to discuss the effect of Higgs 
boson exchange more quantitatively. We take the universal soft masses for the squarks 
and sleptons m , the gaugino mass My2, and the A-terms a at the GUT scale. In order 
to realize a relatively light Higgs mass spectrum, we take the quadratic mass parameters 
of the Higgs potential m 2 H1 and m 2 H2 to be independent. The values of the /i and B 
parameters in low energy are determined to reproduce correct vacuum expectation values. 
Between the GUT and the right-handed neutrino mass scales, the parameters evolve with 
the RGE of minimal SUSY standard model (MSSM) with the right-handed neutrinos, that 
induces (Am- )y which is approximately given in Eq.(^). The evolution of the parameters 

<>L/ 

obeys the MSSM running below the right-handed neutrino scale. In the actual calculation 
we solve relevant RGE's numerically. Once we obtain the low-energy parameters, we 
calculate the coefficients A's and <fs given by the one-loop diagrams [Uj and the Higgs 
exchange effects in Eq.(JHJ). Then we used the the values in Ref. [21] to calculate the 
conversion branching ratios from the coefficients. 

The results of the calculation are shown in FigsH] and 121 We show in FigJIJa) the 
scatter plot of the value of B (// Al — > eAl). We fixed the m and Mi/ 2 parameters to be 
2000 GeV, and tan (3 = 60. The right-handed neutrino masses are also fixed to be M l N = 
10 14 GeV for all the generations, and the Yukawa coupling constants /„ are determined 
so as to fit the neutrino oscillation data of Am^ tm = 3 x 10~ 3 eV 2 , Am 2 ol = 4 x 10~ 5 
eV 2 , sin 2 6» atm = 0.5, and sin 2 9 soi = 0.25 [23 123 1211 • We took tlie normal hierarchy with 
m v \ = and f/ e3 = 0, where U is the Maki-Nakagawa-Sakata matrix [2H1- As expected, 
we can see the enhancement in the light H° region (m#o < 600 GeV) because of the 




(14) 



whereas the contribution from the photon exchange is calculated to be 




(15) 
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Figure 1: The m H o dependences of the branching ratios of the following processes are 
shown: (a) jj, — e conversion in aluminum nucleus and (b) fi —> ey decay. We take the 
right-handed neutrino masses to be 10 14 GeV, and tan/5 = 60. The soft masses for the 
Higgs fields are treated as free parameters. 

Higgs-exchange effect. We can obtain the branching ratios for different values of Mn by 
the scaling of B(fiA\ — ► eAl) oc (Mat) 2 . Fig^b) is the same plot for B(fi —>■ cy). The 
enhancement in the light H° region is absent in the fi —>■ decay. Therefore, taking 
the ratio of the branching ratios (B(fiAl — > eA\)/B(fi — > ej)), we would see whether the 
contribution of the Higgs boson exchange is large. FigOJa) shows this ratio for the cases 
of mo = M1/2 — 500, 1000, and 2000 GeV. The ratio becomes large and can reach 0(1) 
for small m#o. As increasing m#o, it monotonically approaches to 0.0026, which is the 
value predicted for the case that the \i — e conversion occurs through photon exchange 
diagrams j21]. The deviation from 0.0026 indicates that the existence of the operators 
besides the photonic dipole operator. Although this is an interesting prediction of the 
Higgs-exchange LFV, it is not clear whether the Higgs-exchange effect is responsible for 
the deviation when it is measured. 

The target atom dependence of the conversion branching ratio is another interesting 
quantity which can discriminate the different operators of the \x — e conversion process 
|24j . There are three types of operators: dipole, scalar four-fermion, and vector four- 
fermion. For light nuclei, all those operators show similar Z dependences in the \x — e 
conversion amplitude. Namely, the overlap integrals behave like i.e. D/(8e) = = 
yW ~ ((A-Z)/A)S {n) = {{A-Z)/A)V {n \ in the non-relativistic limit of the muon wave 
function. On the other hand, these overlap integrals take different values in the heavy 
nuclei due to the large relativistic effect. For examples, £>(/iPb — > ePb)/B(/j,A\ — > eAl) 
for the dipole, the scalar, and the vector operators are 1.1, 0.70, and 1.4, respectively 
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[21] . Thus the ratio of the branching ratio in a heavy nucleus to that in light one provides 
information on the type of operators responsible for the fi — e conversion. In the case 
where the Higgs-exchange is dominated as in the light H° region, the heavy to light 
ratio corresponds to a value of the scalar-operator prediction, which would be a robust 
indication of the effect. We plot the ratio of _B(/iPb — > ePb)/B(fiAl — > eAl) in Fig|2fb). 
We can see that the values indeed approach to the prediction of the scalar operator (0.70) 
in the light H° region, and they increase and approach asymptotically to the prediction 
of the dipole operator (1.1) as increasing m#o. Measurement of this dependence can 
provide us with useful information which will eventually lead us to extract the size of the 
Higgs-mediated LFV. 

The enhancement due to the Higgs boson exchange is not as significant in \i — > 3e 
process as in /i — e conversion due to the small Yukawa coupling of an electron. The 
ratio B(fi — > 3e)/B(fi — > ej) is shown in Fig|2fc). There we can see the ratio is almost 
constant (0.006) over the same parameter region as in FigJ2fa). 

In summary, we calculated the branching ratios of the coherent fi—e conversion process 
in the SUSY seesaw model. The effect of the Higgs-mediated LFV may give dominant 
contribution in the large tan j3 and the light H° region because of the tan 6 j3 and (m H a)~ 4 
enhancement of the branching ratio. The Higgs-exchange effect gives interesting signals 
on ratios of the branching ratios. For example, the ratio of B(fiAl — > eAl) /B (fx — > ej) 
is predicted to be 0.0026 with only the dipole operator, but the Higgs-mediated process 
enhances this ratio significantly. Moreover, information on the operator responsible for 
the fi — e conversion can be obtained by the target atom dependence of the fi — e conversion 
branching ratio. This is particularly useful because we can make a definite theoretical 
prediction for the ratio like B(fiPh — > ePh) / B(fiAl — > eAl) depending on the dominant 
type of operators. For the ratio of B(fi Al — > eAl) /B (ft — > ej), such definitive prediction 
is possible only in the case of the photon-dipole-operator dominance. 
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Figure 2: The following ratios of the branching ratios are shown as functions of m H o: 
(a) B(fiA\ -> ek\)/B{n -> e7), (b) 5(/iPb -> ePb)/B(//Al -> eAl), and (c) £(// -> 
3e)/B(fi — > e7). We take the right-handed neutrino masses to be 10 14 GeV, and tan/3 = 
60. The soft masses for the Higgs fields are treated as free parameters. 
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